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In numerical evolutions of binary black holes (BBH) it is desirable to easily control the orbital 
eccentricity of the BBH, and the number of orbits completed by the binary through merger. This 
paper presents fitting formulae that allow to choose initial-data parameters for generic precessing 
BBH resulting in an orbital eccentricity ~ 10~*, and that allow to predict the number of orbits 
to merger. We further demonstrate how these fits can be used to choose initial-data parameters 
of desired non-zero eccentricity. For both usage scenarios, no costly exploratory BBH evolutions 
are necessary, but both usage scenarios retain the freedom to refine the fitted parameters further 
based on the results of BBH evolutions. The presented fitting formulas are based on 729 BBH 
configurations which are iteratively reduced to eccentricity 10~*, covering mass-ratios between 
1 and 8 and spin-magnitude up to 0.5. 101 of these configurations are evolved through the BBH 
inspiral phase. 
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I. INTRODUCTION 

The next generation of gravitational-wave detectors, 
such as advanced LIGO, VIRGO and KAGRA [iHl is 
under construction. These advanced detectors will have 
an order of magnitude increase in their sensitivity and 
are widely expected to make the first direct detections 
of gravitational waves by sources such as coalescing com- 
pact object binaries. The ability to detect and study 
these systems depends on the quality and accuracy of 
the theoretical waveform models used in building search 
templates for these detectors. 

After breakthroughs in numerical relativity Q, direct 
numerical calculation of the signal emitted during the last 
phase of the binary's life have become available and have 
been rapidly improving (e.g. 0, Q)- These numerical 
waveforms guide analytical modelers and data analysts 
in the construction of analytical templates [ol-flTj and arc 
used to assess the properties of gravitational wave data- 
analysis pipelines [171 Il8l | . 

The parameter space for numerical studies of binary 
black holes is seven dimensional: Mass-ratio, and the 
two spin-vectors of the two black holes. The total mass 
scales out due to the scale invariance of the vacuum Ein- 
stein equations, and eccentricity is expected to be radi- 
ated awa.y during the preceding gravitational- wave driven 
inspiral jl^, [2^. Moreover, simulations for gravitational 
wave data-analysis need to cover a large number of in- 
spiral orbits and have sufficient accuracy (e.g. the recent 
review [2l|). These requirements increase the computa- 
tional cost of BBH simulations and limits the number of 
simulations that can be performed. Therefore, BBH sim- 
ulations have generally focused on lower-dimensional sub- 
spaces of the entire BBH parameter space, for instance 
non-spinning systems, or systems with spins aligned with 
the angular momentum. For example, the number of dis- 
tinct BBH parameters configurations used in Ninja-2 
was 29: 6 configurations with non-spinning black holes, 
and 20 with aligned spins (some configurations were com- 



puted independently by several groups for a total of 40 
numerical simulations). Numerical simulations of non- 
precessing BBH systems are also much more carefully 
studied, e.g. |2a - l36j . 

Furthermore, most of the work done to compute ana- 
lytical waveform models focuses on non-precessing bina- 
ries [13, [H, [H [IMl- Significantly less work has been 
done in trying to understand dynamics of precessing bi- 
naries and the resulting waveforms, and how the numer- 
ical data agrees with the analytical approximations (45l - 

The first step for evolving a precessing system requires 
low-eccentricity initial data, since the orbit of a isolated 
binary circularizes during the inspiral via the emission of 
gravitational waves [l^ ll^] . Even binaries starting with 
some eccentricity at the beginning of their evolution are 
expected to have a negligible eccentricity near the merger 
phase, when the emitted signal enters the frequency band 
of these ground-based detectors. 

In a numerical relativity simulation of the inspiral of 
binary black holes, the eccentricity is predetermined by 
three free parameters in the initial data: the orbital fre- 
quency r^O: the separation between the holes tq and the 
radial velocity (we often use the dimensionless expan- 
sion factor oo = ro/ro). When initial data is constructed 
with the assumption of circular orbits, the resulting tra- 
jectories have an orbital eccentricity of the order one per- 
cent (53l - [55| . This eccentricity arises by neglecting the 
small, but non-zero, radial velocity and the initial relax- 
ation of the black holes. 

A variety of eccentricity definitions are given in the 
analytical literature |56l-[63| as well as in numerical rel- 
ativity [ssl - lssi I63l466l |. All of these definitions employ 
residual oscillations in the orbital variables such as the 
orbital frequency, proper horizon separation and coor- 
dinate separation to estimate the eccentricity. Several 
methods are introduced to choose initial-data parame- 
ters that result in lower eccentricity. The evolution of 
post-Newtonian equations was used to find quasi-circular 
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parameters for the trajectories of binaries 63| . More re- 
cently, iterative procedures were developed to remove ec- 
centricity from the initial data [13, EB-E^ ■ These meth- 
ods first utilize post-Newtonian (PN) information to find 
initial data with reasonably low eccentricity. The ini- 
tial data is evolved for about two to three orbits, and 
after analyzing the orbit, the initial-data parameters arc 
corrected to reduce the orbital eccentricity. Using these 
new initial-data parameters, the procedure is repeated 
until the desired value of eccentricity is obtained. For 
non-precessing binaries, the method using proper hori- 
zon separation worked well in reducing the eccentricity 
of binaries [13, [H, [t^. For precessing binaries, a new 
method was introduced in Ref. which uses the in- 
stantaneous orbital frequency to reduce the eccentricity 
to below lO"*^ in about four iterations. 



Iterative eccentricity removal works well j6a |. however, 
it introduces extra steps into the numerical simulation of 
BBH systems. The goal of this paper is to deal with ec- 
centricity removal for conformally flat BBH systems once 
and for all: We systematically apply the eccentricity re- 
moval procedure to a large number of different BBH sys- 
tems of different mass-ratios, spin-magnitudes and spin- 
directions. We cover mass-ratios between 1 and 8, and 
spin- magnitudes up to 0.5. We then perform fits that will 
allow us to predict low-eccentricity initial-data parameter 
for BBH systems with any spin-orientations. As we show, 
these fits result in initial conditions with remaining ec- 
centricity of ~ 10"''. This is quite likely sufficient for all 
near-term GW data-analysis purposes 
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eccentricity fits can also be used to predict initial data 
parameters of a desired non-zero eccentricity. Finally, 
we evolve 101 of the 729 low-eccentricity initial data sets 
through their entire inspiral, and use this information to 
prepare a fitting formula that allows us to predict the 
number of inspiral orbits for generic BBH binaries. 



This paper is organized as follows: In Sec. |TT1 we dis- 
cuss eccentricity and motivate a new definition of eccen- 
tricity, cr, that uses the radial velocity in addition to 
the orbital frequency. Once low-eccentricity parameters 
are known, allows us to estimate the eccentricity of 
any parameter choice nearby without evolving the initial 
data. In Sec. lIIIl we quasi-circularize iteratively 729 BBH 
configuration with different mass-ratios, spin-magnitudes 
and -orientations, and various initial separations. Three 
of these simulations are used to compare our newly de- 
fined eccentricity cr to the standard eccentricity based on 
the orbital frequency, en, used in Ref. [i^. Section ITVl 
introduces fitting formulae to the low-eccentricity con- 
figurations, and demonstrates their efficiency and their 
advantages over post- Newtonian formulas. In Sec. El we 
provide formulas for generating initial-data parameters 
that result in evolutions of a predetermined eccentric- 
ity, and assess the quality of these formulas. Finally, we 
summarize our conclusions in Sec. IVIl 



II. DEFINING AND REDUCING 
ECCENTRICITY 

A. Newtonian dynamics 

For two bodies on a Newtonian orbit with small eccen- 
tricity e, the distance r(t) between their centers and the 
orbital frequency fl{t) can be written as 

r{t) = f [1 -e sin(17t4-0o)] +C(e^), (1) 
n{t) = Cl [l + 2esm{nt + (l3o)] +0{e'^). (2) 

Here r is the average separation, D, is the average orbital 
frequency and (j)Q is a phase component. Taking a time- 
derivative of Eqs. [T] and [2l we find 



r{t) 



-feil cos(m + (j>o) + 0(6^), 
2en^cos{nt + (j)o) + 0{e^). 



(3) 
(4) 



Note that the factor f2^ is the product of the average or- 
bital frequency Q and the frequency of the oscillations in 
fl{t). In Newtonian gravity, these two frequencies agree. 

For circular orbits, the distance between the Newto- 
nian masses tq and their orbital frequency are related 
by Kepler's law. 
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where m is the total mass of the binary, and G represents 
Newton's constant which we will henceforth set equal to 
unity. In this case, eccentricity e and radial velocity r{t) 
are both identically zero. 

Let us consider how small perturbations in orbital sep- 
aration, orbital frequency and radial velocity will affect 
the eccentricity of the originally circular orbit. To accom- 
plish this, we relate a circular orbit with tq and JIq with 
a slightly eccentric orbit given by Eqs. ([Ij and ©. At 
time f = 0, we set r(0) = ro -I- At, il{0) = ^0-1- Afi, and 
wc set a radial velocity f(0) = Af. The perturbed orbit 
will in general have an average distance f different from 
To, and an average orbital frequency different from flo- 
Wc write 



ro + Sf, n^Qo + <5n. 



(6) 



Kepler's third law Eq. ([5]) is also valid for r and fl, and 
from it follows 
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The change Sf and SQ are thus determined by one small 
parameter e, which is as of yet undetermined. To pro- 
ceed we substitute Eq. © into Eqs. (H])-®, evaluate 
these expressions at i = 0, and equate to the assumed 
perturbations: 

ro + Ar = (l + e)ro(l-esin0o)'^o, (8) 
Qo + An = (l-3e/2)(l + 2esin0o)f^o, (9) 
Af = -(l + e)ro(l-3/2e)17oecos0o- (10) 
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To first order in e and e, these equations simplify to 
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-j'o ^^0 e cos 00- 
Inverting Eqs. ([TTll-dlSl) yields 
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To summarize, perturbing a circular Newtonian orbit by 
Ar, Afi and Ar results in an orbit with eccentricity 
the perturbed orbit has an average radius larger by 
Sf/r = e relative to the the circular, unperturbed orbit. 



B. Numerical relativity 

Let us now consider general relativistic BBH bina- 
ries computed by numerical simulations. As in earlier 
workjsBi l66l| we define eccentricity based on periodic os- 
cillations of the separation of the binary or the orbital fre- 
quency, carrying over the Newtonian definitions. Specif- 
ically, the instantaneous orbital frequency J7(t) and the 
distance r{t) are calculated from the coordinate motion 
of the apparent horizons' centers. We define their rela- 
tive separation vector r{t) = Ci{t) — C2(t) with magni- 
tude r{t) = \r{t)\, where Ci{t) are the coordinates of the 
center of each black hole. Using standard Euclidean vec- 
tor calculus, the instantaneous orbital frequency is then 
computed as 



rjt) X m 

r^{t) ' 



(17) 



where Cl(t) is its magnitude. 

A compact binary inspiral starts at i = with an ini- 
tial separation ro, an orbital frequency and a radial 
velocity fo. The time derivative of the orbital frequency 
Cl{t) is computed and fitted with the functional form 
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The first part of the fit, S'f2,fit(i), is a non-oscillatory func- 
tion that captures the radiation-reaction driven inspiral. 
The second part captures the oscillatory contribution of 
the orbital eccentricity using fitting parameters i3a, 
and (j). Note that the sinusoidal frequency in Eq. [TH] is 
given by the radial frequency fl^. In general relativity, 
the orbital frequency exceeds the radial frequency fir 



causing periastron advance [63. With the fitted param- 
eter equal to 2eflflr (cf . the comment after Eq. |4]) , 
we can derive updating formulas based on 0(t) [66l |: 
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where Bq and (j> are fitted for while J7o and ro are given 
by the initial data. In Eq. we have introduced the 

expansion factor clq = ro/rQ which appears naturally in 
our formulation of the initial-value problem for BBH bi- 
naries with radial velocity [s^ . We use Eqs. (|19p and ((20)) 
to reduce eccentricity iteratively for every simulation in 
this paper. Using the fitted Bq and fir in addition to 
r^o, we estimate the eccentricity en as 
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Equation ([14)) gives the eccentricity of a Newtonian or- 
bit that differs by Aft and Ar from being circular. If we 
know non-eccentric initial data parameters D,o.e=o and 
flo,e=0: we can use the deviations from the e param- 
eters to define eccentricity: 



0-0 ~ ao,e=o 



^0,0 — flo,e=Q 



1/2 



(22) 



We will test this formula in Sec. lIIIl and use it in Sec.fVlto 
propose a technique to construct BBH initial-data with 
specified non-zero eccentricity. 



III. LOW-ECCENTRICITY BBH PARAMETERS 

Our goal is to compute several hundred low eccentric- 
ity initial-data sets for subsequent evolutions. We pro- 
ceed in two stages: In stage 1, we circularize 69 binaries 
by starting with post-Newtonian estimates for the ini- 
tial orbital frequency and radial velocity. We perform 
three to four levels of iterative eccentricity removal. We 
use the low-eccentricity initial data parameters to gen- 
erate fitting formulae that predict low-eccentricity initial 
data parameters for other choices of masses and spins. 
In stage 2, we use these fitting formulae to compute 660 
additional BBH configurations, which we find have an 
average eccentricity of about 0.003, i.e. a factor of ~ 5 
below the first stage simulations. We finally apply itera- 
tive eccentricity removal to the stage 2 simulations until 
the eccentricity is less than 0.001. 



A. Numerical methods 

We construct BH-BH initial data using the conformal 
thin sandwich formalism (72l . pTSj with quasi-equilibrium 
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FIG. 1: Eccentricity removal based on time derivative 
of the orbital frequency dil/dt, applied to a precessing 
binary black hole, ScN 21. Shown is Cl vs. time and four 
eccentricity-removal iterations. 



boundary conditions j74l - l76l | and a radial velocity as in 
Ref. [13] ■ AH initial-data sets considered here use confor- 
mal flatness and maximal slicing. The resulting set of five 
nonlinear coupled elliptic equations is solved with multi- 
domain pseudo-spectral techniques described in Ref. ITTII 
implemented in the Spectral Einstein Code SpEC [78| . 
Calculation of initial data with desired physical param- 
eters (masses mA,?^^s and dimensionless spins XA, Xb) 
requires a root-finding procedure to determine the cor- 
responding initial data parameters (radii of the excision 
boundaries ta, rs and angular velocities of the horizons, 
a A and fJs). This root-finding is described in Ref. [s^. 

The initial data is evolved with SpEC [tI], using a 
first-order representation (79j of the generalized harmonic 
svstem 16 . ISOl l8ll | that includes constraints damping 
terms p, l79l l82l . Constraint damping parameters are 
chosen as in [36| . based on experience gathered in (30| . 
The computational domain extends from excision bound- 
aries located just inside each apparent horizon to some 
large radius, where the outgoing gravitational radiation 
pass freely through the outer boundary. Outer boundary 
conditions [t^, HH, [s^I are iinposed to prevent the influx 
of constraint violations and undesired incoming 

gravitational radiation j92l l93j| . while no boundary condi- 
tions are imposed at the inner excision boundaries. Inter- 
domain boundary conditions are enforced with a penalty 
method [§3, [sEl • The overall evolution techniques (con- 
straint damping parameters, choice of domain dccompo- 
sition) are essentially identical to the inspiral phase of 
Ref. [36|. For precessing runs, we employ a coordinate 
transformation based on pitch- and yaw-angles (i.e. two 
of the Euler angles) , described in detail in a forthcoming 
publication [o^ . This technique works well for moderate 
precession as in the cases presented here, although in fu- 
ture runs, it will be replaced by the more sophisticated 
coordinate transformations developed in (96j . 
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FIG. 2: Convergence of the eccentricity-removal pro- 
cedures in the {SQ,e, Sae) plane. In the upper panel, we 
show the eccentricity removal sequence of Fig. [1] The red 
number next to each point gives the eccentricity en, while 
the blue number next to each blue solid line gives the relative 
eccentricity eR. In the lower panel, we plot the eccentricities 
ratio en/eR. 



B. Iterative Eccentricity Reduction 

We shall start by describing eccentricity removal for 
one typical precessing configuration, labeled ScN 21. 
This configuration has a mass-ratio of g = 771^1/7713 = 1.5 
and both black holes have dimensionless spins 0.5 ini- 
tially tangent to the orbital plane. The spin of the 
larger black hole points exactly away from the smaller 
black hole, whereas the spin of the smaller black hole 
is anti-aligned with its initial velocity. The initial co- 
ordinate separation between the holes is r = 16m. We 
begin eccentricity removal with orbital parameters de- 
termined from a non-spinning PN approximant to choose 
r^o 0.014427/771, and we use a previous good approxi- 
mation for do = —3.6 X 10~^ from other simulations. 

Wc evolve the binary with the initial orbital parame- 
ters {Qo,do) for about two orbits, and record the time- 
derivative of the orbital frequency, Q(t). This data is 
plotted in Fig. [1] as "Iter 0" and we deduce an eccen- 
tricity of en ~ 0.018. Equations ^ and ^ give the 
improved values for Qq and do to use in next iteration, 
labeled "Iter 1" in Fig. [T] with eccentricity en = 0.0037. 
This procedure is repeated three more times until a final 
eccentricity en ~ 1.5 x 10"'* is achieved. Note that as 
the eccentricity falls below e ^ 0.001, spin-induced oscil- 
lations become a ppa rent at twice the orbital frequency, as 
discussed in Ref.[66[. The spin-induced oscillation domi- 
nate for e ^ 0.0003, making unambiguous determination 
of a residual orbital eccentricity difficult. 

Equation indicates that the eccentricity is the 
square sum of two components, with one measuring the 
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FIG. 3: Distribution of configurations in set iSo. Given 
are the number of non-spinning (black), spinning but non- 
precessing (red) and precessing binaries (blue) for different 
mass-ratio q. 



needed change in orbital frequency and the other the 
change in radial velocity (or derivative of expansion fac- 
tor). When plotting parameters during eccentricity re- 
moval, Eq. ((22|) suggests a natural scaling of the axes, 
relative to the configuration with lowest eccentricity: 
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In these axes, the Euclidean distance to the origin should 
correspond directly to the eccentricity, cf. Eq. (fT4)) . Fig- 
ure [2] shows the sequence of the eccentricity reduction 
iterations shown in Fig. [T] plotted in these coordinates. 
The upper panel of the plot shows SQ^ and dd^ as the 
binary is quasicircularized itcratively. The distance in 
the ((5de,5rie) plane between any point and the origin 
corresponding to "Iter 4" is equal to the eccentricity e^ . 
Because Iter 4 is not exactly at zero eccentricity, this 
distance is only an approximation to Eq. (|14p . thus ex- 
plaining the raise of en/ en at iterations 2 and 3 in the 
lower panel of Fig. [21 



C. Binaries in data set So 

As a first step, we quasi-circularize itcratively a set 
of 69 different configurations, which we shall refer to as 
set Sq. The binaries in iSg have mass-ratios 1 < < 8. 
There are 9 non-spinning binaries, 24 spinning binaries 
with aligned spins (i.e. without precession), and 36 pre- 
cessing binaries. The dimensionless black hole spin is 
generally 0.5, but sometimes smaller. The parameters 
for all 69 configurations are given the first nine columns 
of Table U in the Appendix, and Fig. [3] show the distri- 
bution of non-spinning/ aligned spin/non-aligned spins as 
a function of mass ratio: Most precessing binaries have 
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FIG. 4: Convergence of eccentricity estimates en and 
er for cases ScN 21, 50 and 68. Both eccentricity defi- 
nitions decrease at the same rate to a value below the 10""^ 
limit. The parameters of the three configurations are listed 
in Table [J 



mass ratio either 1.5 or 3. At mass ratio 8, no precessing 
binaries are evolved. 

Wc apply the iterative eccentricity reduction method 
to all binaries of set Sq. Figure |4] illustrates performance 
of the iterative eccentricity removal for three cases. This 
figure compares also the eccentricities eo and cr. Both 
eccentricities decrease at the same rate, but the estimated 
Cfj is always smaller than eji, since it is estimated rela- 
tive to an iteration with finite, albeit small, eccentric- 
ity. As we discuss below, is a useful measure to es- 
timate the eccentricity of runs with the same physical 
parameters but differing orbital parameters flo and do, 
once the corresponding quasicircularized orbital param- 
eters r2o,e=o and do,e=o are known. 

Application of iterative eccentricity removal to all bi- 
naries in So results in orbital parameters VIq and do as 
listed in Tabic HI with an estimated orbital eccentricity 
CQ. Eccentricity removal is terminated once an eccen- 
tricity eo < 10"'^ is reached except for few runs that 
are intentionally left with a larger eccentricity for evolu- 
tions to compute periastron advance (cf. [40j). The initial 
and final eccentricities of binaries in set So are shown in 
Fig. [51 The red points show the initial eccentricity while 
the blue point correspond to the final eccentricity. 

The initial eccentricities for most of these runs are 
larger than 0.01. The initial conditions for these runs 
used FN approximants TaylorT3 to predict the values 
of the orbital frequency at a given separation, and the 
radial scaled velocity is set to a small constant value 
(~ —3 X 10~^). However, for a few runs, the initial ec- 
centricities are much smaller, almost 10"'^. For these 
simulations, we used Kepler's law to predict the initial 
orbital frequency given the quasi-circular parameters of 
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FIG. 5: Initial and final eccentricities of the configu- 
rations in set So- The eccentricity reduction target was 
10^^, except for runs intentionally left at larger eccentricity 
for calculations of periastron-advance. 



a closely similar configuration. The procedure based on 
Kepler's law was experimental and motivated the devel- 
opment of more rigorous fitting formulas to predict the 
initial parameters of quasi-circular binaries, as presented 
in later parts of this paper. 

All quasi-circular binaries in Sq were evolved through 
the inspiral phase until approximately 1 — 2 orbits before 
merger. To provide the reader some context about the 
length of these simulations, Table U lists also the number 
of evolved orbits up to a orbital frequency Mfi = 0.05. 
For informational purposes, we also report in TableUthe 
orbital frequency mil f at which the simple inspiral evolu- 
tion techniques failed and the number of orbits Nf com- 
pleted by then. Seven binaries were evolved for more than 
30 orbits. The evolutions of merger and ringdown of these 
binaries require significant extension and refinements of 
previous spectral merger-techniques [s^ |93] , and will be 
reported elsewhere. 



D. Binary configurations Si, . . . , S5 

As Fig. [S] demonstrates, iterative eccentricity removal 
by brute force is possible. However, even for some config- 
urations of Fig.[SJ the use of ad hoc fitting formulas based 
on Kepler's law resulted in significantly improved initial 
guesses fig and do . This allowed to begin iterative eccen- 
tricity removal at smaller initial eccentricity, reducing the 
number of eccentricity-removal iterations and thus low- 
ering the computational cost of iterative eccentricity re- 
moval. The goal of this section, therefore, is to construct 
a much larger sample of low-eccentricity orbital configu- 
rations, which will be used in Sec. IV to develop fitting 
formulae valid for generic BBH systems. 

In addition to the 69 configurations of 5o, we circular- 
ized a total of further 660 configurations, of which 156 
are non-precessing (of which 40 are non-spinning) and 
504 are precessing. The configurations fall into five sets. 
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FIG. 6: Parameter space covered with eccentricity removal. 
Each point in set <Si represents one configuration for a total of 
40. Each point in set S2 represents five different xa directions, 
for a total of 190. Each point in set S3 represents 15 different 
XA and xb directions, for a total of 300. The set 54 and 
55 contain 130 configurations with randomly generated spin 
directions and (for 32 of these) randomly generated q. 



as indicated in Fig. [HI 

• Set Si contains 40 non-spinning BBH configura- 
tions. 

• Set 1S2 contains 190 single-spin configurations 
where the more massive black hole carries spin 
Xa = 0.5. For each {q^r/m) pair, five dif- 
ferent spin-directions are considered, with 9a = 
0, 7r/4, 7r/2, 37r/4, tt. The first and the last of these 
9a represent spins aligned and anti-aligned with the 
orbital angular momentum, so that 76 configura- 
tions are non-precessing, while the remaining 114 
configurations are precessing. 

• Set iSs contains 300 configurations where both 
black holes carry spin XA = Xb = 0.5. We consider 
the same five spin directions for the more massive 
black hole as in set ^2. For each of these, we con- 
sider three spin directions for the less massive black 
hole: 9b = 0,7r/4,7r/2. 40 of the configurations in 
set S3 are non-precessing, and the remaining 260 
are precessing. 

• Set ^4 contains 98 binaries with spin-magnitudes 
XA = Xb = 0.5, and random spin-directions. 25 bi- 
naries were selected for each combination of mass- 
ratio q — 1 or q ~ 2, and of initial separation 
r/m = 16 or r/m = 17. Two of these runs were 
compromised, so that this results in 98 configura- 
tions. 

• Set 1S5 contains 32 binaries with random spin- 
magnitudes XA,B < 0.5 at random spin-directions 
with random mass-ratios g € [1,2] for fixed r/m = 
15. The explicit parameters are listed in Table [H] 
in the Appendix. 
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FIG. 7: Initial and final eccentricities of all runs in 
the paper. The first 69 configurations corresponds to tlie 
set So of binaries listed in Table [J The next 530 configura- 
tions corresponds to sets Si,S2,S3 of binaries, and the last 
130 configurations correspond to sets S4, and S5 with random 
spin configurations. The final eccentricity of the runs in sets 
Si, . . .St, is reduced to less than 10^^. 

Figure [7] summarizes initial and final eccentricity of all 
configurations which were quasi-circularized. We initial- 
ized eccentricity reduction using an ad-hoc formula based 
on Kepler's law. As can be seen in Fig. [71 this reduced the 
starting eccentricity of 51^2.3 by about a factor of 10 rela- 
tive to that of 5o. However, it did not significantly reduce 
the starting eccentricity of the processing systems in 54,5. 
A more comprehensive approach is therefore needed. 



IV. ECCENTRICITY REDUCTION USING 
FITTING FORMULAS 

Kepler's law helped in quasicircularizing the binaries 
of the set So when small variations are introduced in 
either the separation or the mass ratio of the binary. The 
Newtonian Kepler's law does not incorporate spin effects 
in the initial parameters of the binary. Therefore, we 
will base fitting formulas for spinning binaries on post- 
Newtonian expansions for spinning binaries. 



A. PN formulas 

In the post-Newtonian approximations, the initial 
quasi-circular parameters fig or '^o (and rg) are given 
as a series expansion in m/r (or in fio) for any physical 
configuration for the binary. We shall use the PN expres- 
sions derived by Kidder |98| . In these expressions, non- 
spinning effects are included up to 2nd post-Newtonian 
order, while the spin effects enter at 1.5th and 2nd 
post-Newtonian order. We define the unit vector along 
the angular momentum Ln = Ln/|Ln| where Ln = 
/i(r X v) and the spin unit vectors sa,b — ^a/b/ Sa/b = 
Xa/b/xa/b- Furthermore, let m = Ma + Mb be the 
total mass, q = mA/iriB the mass-ratio, = mAmB/m 
be the reduced mass and ij = mAmB/rm? = q/{l + g)^ 
the symmetric mass-ratio. In terms of these quantities, 
Kidder derives the orbital frequency as a function of 
separation. 
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The radial velocity as a function of separation is given by: 
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The orbital separation r/m as a function of orbital frequency is 
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And finally, the number of orbits accumulated during inspiral from initial orbital frequency Vli to final frequency Q / 
is given by = [^'(rii) — 'I'(ri/)] /{2tt) with orbital phase ^' satisfying 
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B. NR Fitting Formulas 

The goal of this section is to provide fitting formulae 
that can predict f^o and do as functions of symmetric 
mass-ratio rj, initial black hole spins Sj^/b and initial 
black hole distance r/m. This is an 8-dimensional pa- 
rameter space. As we shall see, linear fits are not suf- 
ficient, however, if one were to write a straightforward 
8-dimensional higher-order Taylor-series, the number of 
coefficients would quickly explore with expansion order. 
Therefore, one needs to be judicious to include only terms 
that contribute to the fit. 

Motivated by Eqs. ([25]) and ((26)) . we do not write down 
fitting formulae directly for flo and do, but rather for the 
quantities 



n2 3 

"o'o 



The functions 6ks and 5pg should be functions of the 
spins of the black holes, mass-ratio and separation. These 
functions should vanish for zero spins, to reduce to the 
fit for non-spinning BBH. 

We will use in total 14 basis-functions to represent the 
spin-sector, which we shall label Ba, a = 1, . . . , 14. The 
first four basis-functions take the functional form of the 
spin-terms in Eqs. (j25p 



(35a) 
(35b) 
(35c) 
(35d) 



(=1 = Sa • L, 
62 = Sb • L, 

(34 = Sa,± ■ Ss I , 



(29) 



and 



Po 



5m /ro 
6477 \m 



ao . 



(30) 



For each quasi-circular configuration, we compute ko and 
Po- 

Beginning to fit the non-spinning components, we only 
consider the quasi-circular configurations of set iSi and fit 
polynomials in t] and inverse distance u = m/r: 



i+j<3 
i+j<3 

pm = 



(31) 
(32) 



(Here, the subscript 'NS' refers to non-spinning bina- 
ries). We found that the triangular truncation i -I- j < 3 
gave a good fit with a reasonably low number of coeffi- 
cients. The coefficients a,; j and bij are listed in Tables Hill 
andHVl 

For spinning binaries of sets ^2, S3, S4 and S^, we first 
subtract the fits Eqs. ([51]) and ([5^ . i.e. we compute for 
each spinnning quasi-circular configuration 



Sks = Ko — Kjvs|r;o,«o ' 

Sps = Po - PNslm.uo ■ 



(33) 
(34) 



where Si^± = Si — L(L • S^), i = A,B are the projections 
of each black hole spin into the orbital plane. 

The next two basis-functions account for the the angle 
between Sa/b relative to the other black hole: 



65 

66 



(35e) 
(35f) 



where f is the unit-vector pointing from black hole 2 to 
black hole 1. These terms are not present in the post- 
Newtonian expansion. However, we find them to be nec- 
essary for a good fit. We conjecture that the origin of 
these terms rests in details of how numerical BBH initial 
data is constructed. For instance, conformally flat quasi- 
equilibrium BBH initial data requires a somewhat dif- 
ferent choice of do to achieve low ecc entri city compared 
to superposed Kerr-Schild data llOd ]. Both types 
of initial-data require some offset in do relative to post- 
Newtonian expansions, and the need to include £=5 and eg 
when fitting for low-eccentricity initial data parameters 
indicates that this offset depends on the angle between 
black hole spin and direction to the partner black hole. 

Finally, as demonstrated below, the fit of the spin- 
sector can be improved by including certain higher-order 
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FIG. 8: Eccentricities en and er for the binaries of tiie sets 
52, and those configurations of 53,4,5 with en < 0.0001. The 
residuals er of tlie fitting formulas are comparable to the 
actual estimated eccentricity en. Further, improvement of 
the fitting formulas requires quasi-circular configuration with 
lower eccentricities en. The first 40 simulations are nonspin- 
ning. The expected PN eccentricities are plotted in grey. 



basis- functions: 



es 
eg 
eio 
en 
ei2 
ei3 
ei4 



(Sa-L)2, 
(Sb-L)2, 

(Sa 
(Sa, 

(Sb 



SB)^ 



(35g) 
(35h) 
(351) 
(35j) 
(35k) 
(351) 
(35m) 
(35n) 



The functions Sks and Sps arc expanded into a scries 
that depends on 77, u and the spin basis ia- 



^i,j^aV ^ ea , 



i=0 j = 2 q;=1 

3 5 14 

5ps = E E E di,j,artu'''^ec 

i=0 j=2 a=l 



(36) 



(37) 



These expansions have 224 coefficients each. For our best 
result, we fit these coefficients against that subset of con- 
figurations of the sets 53,4,5 which have < 10^*. These 
configurations are located in the green shaded region of 
Fig. \7\ The resulting coefficients arc listed in Tables IVl 
and IVII in the Appendix. 



C. Quality of the fitting formulas 

The quasi-circular orbital parameters that form the ba- 
sis of the fits Eqs. ([55]) and (|37p arc not perfectly accu- 
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FIG. 9: Performance of fit vs. the eccentricity of the 
configurations used in the fit. Eccentricity fits are per- 
formed only using those configurations with en < en, max. As 
a function of en, max the plot shows the average and the median 
of e_R. Reducing the eccentricity en, max improves the fitting 
formulas performance. The numbers indicate the number of 
configurations used in each fit. 



rate; rather they contains residual eccentricity as indi- 
cated in Fig. [T] Therefore, we should expect the low- 
eccentricity fits to be good to only a comparable level of 
precision. If the fits were more accurate than indicated 
by Fig. [3 this would indicate a fit to the errors in the 
numerical data, which would not carry any additional 
physical information. To quantify the quality of the fits, 
we translate the residuals 



Ado = Go, NR. — ao,fit 
into an equivalent eccentricity via Eq. (|15p : 



Ado 



, Ario 

'no. 



NR. 



1/2 



(38) 
(39) 



(40) 



Our target is that cr as computed from the fits is roughly 
comparable to the eccentricity en of the data used to 
compute the fits. If br ^ eo the fit is not as accurate 
as it could be; if ej^ ^ e^, the fit has so many degrees 
of freedom that it fits the residual eccentricity of the NR 
simulations. Figure [5] demonstrates that indeed ~ en, 
as desired. 

The dependence of the quality of the fit onto the con- 
figurations used for the fit is further explored in Fig. |9l 
For this figure, we perform the fit for the spinning sec- 
tor 6ks, 6ps several times, based on those configurations 
within 53,4,5 with eccentricity e^ < en,nia.x- The dashed 
curves of Fig. [5] show the residual eccentricity e^ of these 
fits as a function of threshold eccentricity en. max- For the 
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FIG. 10: Eccentricities {Sae,5^1e) for the binaries of 
the sets Si,-- - ,£5 with e < 0.0001 as predicted by the 
the fitting formulas and the PN expressions. Most binaries 
have an eccentricity less than 0.0001 when using the fitting 
formulas, while PN expressions generate binaries with nearly 
0.01 eccentricity. The fits yield lower eccentricity than PN 
expansions, with the improvement primarily due to the more 
accurate prediction of radial velocity (Sae). We plot the fits 
based on binaries with eccentricity less than 0.0001. 



Because we know low-eccentricity orbital parameters 
^^o,NR, oo.NR these formulae allow us to estimate the ec- 
centricity the PN parameters would have without having 
to evolve with the PN initial-data parameters. Figure [5] 
shows this estimated eccentricity e^ pN as the grey dots. 
Our fits improve the post-Newtonian initial data param- 
eters by about two orders of magnitude, consistently for 
all spin-directions, including precessing binaries. 

Figure [10] plots separately the two components 6ile 
and 6de, cf. Eqs. (jH)) and The plot shows the PN 

data in grey and the results of the fitting formulas as red 
circles. Once again, the improvement gained by the fits is 
apparent: The eccentricity is smaller by nearly two order 
of magnitude relative to PN initial-data parameters. We 
also note that for the fit, Sde and Sfle have approximately 
equal magnitudes. For PN, in contrast, Sfle is about a 
factor 10 larger than ddg, indicating that an error in 
the orbital frequency is the major cause of eccentricity 
when using PN parameters. 



D. Number of orbits fits 



dashed lines the residual of the fit is similar to eo^niax, 
indicating that the fit works as well as can be expected. 
If the basis-functions for the spin-sector are restricted to 
only the linear terms, Cq,, a = 1, . . . , 6, then we obtain 
the dotted lines. Now the residual eccentricity is sev- 
eral times larger than en, max, indicating the possibility 
to improve the fit. 

Wc now compare our low-eccentricity fits with post- 
Newtonian expansions that predict initial data parame- 
ters. Using Eqns. ([^5)) and ([25)1 we compute fipN and 
dpN for each configuration. We then compute the dif- 
ference between eccentricity-reduced NR parameters and 
the PN parameters. 



OO.PN — ao.NR 



sn^ EE 2 

Scie = 

efl,PN = + Sal 



no. 



NR 



(41) 

(42) 
(43) 



When selecting initial data parameters for an evolu- 
tion, it is useful to be able to estimate of how long the 
simulation will be. One widely used measure is the num- 
ber of orbits the binary completes before merger. We 
have evolved the binaries in sets Sq and S5 through in- 
spiral phase, so we can use this information to derive a 
fitting formula for the number of orbits. We base this 
fitting formula on Eqs. (pS)) . specifically 



iV= — (vl/fit(r!,)-vl/g,(r!y)). 



(44) 



The fitting parameters di, . . 
phase-formula. 



. dg are incorporated into the 
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FIG. 11: The number of orbits of the runs and the 
residual of the fitting formulas. In the upper panel, we 
plot the number of orbits A^o.os through the inspiral in red 
between the initial orbital frequency ili and the final orbital 
frequency Q = 0.05/m evolved. In the lower panel, we plot the 
difference between the fitting polynomial and the numerical 
data of the number of orbits in red. We compare the results 
to the PN approximation ploted in orange. 



10" I 
10'" 
C3 10" 



10" 



o 


ScN 1 


o 


ScN 12 



^8 

10"' ^, o,, 



o 



10"^ 



10" 



I5rl/r 



10"^ 



10"^ 



10" I 
10 



10"^ 



^"''ro 



o X- 



10"^ 



10" 



10"^ 

I5a|/a 



10"' 



10"' 



FIG. 12: Eccentricity of runs with initial data parameters per- 
turbed from their low-eccentricity values. Two representative 
configurations are shown, ScN 1 (equal mass, zero spin), and 
ScN 12 (q=1.5, precessing spin). Circles are for positive vari- 
ations, while crosses are for negative variation. The dashed 
line represents the expectation, Eq. (|14[). 



The number of orbits is fitted to Eqs. and (05]) 

using the inspiral data from sets So and ^5 as given in 
Tables HI and im We include only those inspirals that reach 
orbital frequency 0.05/m or greater, that inspiral for at 
least 10 orbits and that have an eccentricity of less than 
0.005. 

The results of the fit for N are shown in Fig. [TTJ The 
upper panel of Fig. [TT]plot the number of orbits evolved 
between the initial orbital frequency fi^ and the final or- 
bital frequency 0.05/to for the inspirals used in these fits. 
For the first 69 simulations of set Sq, binaries with large 
differences in the initial separation and mass ratios are 
evolved. Therefore, a large variation in the number of 
orbits is observed. For the last 32 runs of set 1S5, all bi- 
naries start with the same initial separation but with a 
random mass ratio and spin orientation and magnitude. 
In this case, their number of orbits varies between 15 and 
20. 

In the lower panel of Fig. [Til we plot the difference 
between the number of orbits measured up to a certain 
orbital frequency, and the number predicted by the fit- 
ting formula ((44)) as red points. This residual is less than 
1 orbit for the set Sq, and it is less than 0.2 orbits for 
the binaries of set 55. No variation in the quality of the 
fits is noticed as we changed the final frequency from 
0.05/m. We use also the PN expression in Eq. [28]to plot 
the difference in the number of orbits between the nu- 
merical simulations and the PN predictions measured at 
an orbital frequency of 0.05/to in orange. The initial fre- 
quency used is the value for which the simulation starts 
with, i.e. Hi. The straight post-Newtonian formula pre- 
dicts a number of orbits that generally differs by 3 to 6 



orbits from the NR result. 



V. GENERATING INITIAL DATA WITH A 
PRE-DETERMINED ECCENTRICITY 

So far, we have been concerned with achieving BH-BH 
simulations with very small eccentricity. Let us now con- 
sider how to choose initial-data parameters that result 
in some desired non-zero eccentricity. If one starts from 
known zero-eccentricity orbital parameters (ro, f2oj cto)j 
and perturbs these, (tq + Sr, Qq + SQ, cto + Sd), then the 
Newtonian formula Eq. (|14p should give a reasonable ap- 
proximation of the resulting eccentricity. To test this 
assumption we use two eccentricity-removed configura- 
tions, ScN 1 (equal mass, zero spin) and ScN 12 (mass- 
ratio 1.5, precessing). We perturb first by changing the 
initial separation while keeping Qq and do constant, and, 
second, by changing Qq while keeping dg and do constant. 
For each perturbation, we solve for a new initial-data 
set, and evolve long enough to measure en- The results 
are plotted in Fig. [T^ along with the result of Eq. (|14p . 
e = 3|i5r|/r and e = 2|(5r2|/ri, respectively. 

The agreement between the full numerical simulation 
and the Newtonian formula is very good. The exceptions 
are very low eccentricity en ^ 10~^, where measuring 
eccentricity in the numerical simulation is difficult, and 
for very large eccentricity e ~ 0.1, where linear pertur- 
bations may no longer be adequate. Besides controlling 
eccentricity, one can also use Eq. (|16p to control the phase 
or periastron. 
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So far, we have perturbed off a zero-eccentricity con- 
figuration. As we have seen in Sec. Illli our fitting for- 
mulae for initial data-parameters result in eccentricities 
e ^ 10^^. If one is interested in eccentricities larger than 
this value, then one can also perturb around the results 
of the low-eccentricity fitting formulas. This now allows 
to obtain a BH-BH simulation with desired length (via 
Eq. [^5]) and desired eccentricity, via the fitting formulas 
for low-eccentricity parameters and Eq. (|T4)) and (IT6l) . 
without any iterative procedures. 



VI. CONCLUSIONS 

Evolutions of binary black holes start from initial data. 
Choices that enter the initial data determine the orbital 
eccentricity of the subsequent evolution, and the length 
of the inspiral (i.e. the number of orbits to merger). 
This paper presents techniques that allow to choose ini- 
tial data parameters that result in very small eccentricity 
(e ~ 10~*) with an inspiral of a desired number of or- 
bits. We also present techniques that allow to choose 
initial data parameters that result in a desired non-zero 
eccentricity. Our techniques are applicable for generic 
precessing binary black holes of mass-ratios q ^ & and 
for spin- magnitudes < 0.5. Because of fitting formulae 
we develop here, the use of these techniques requires no 
significant computational cost. 

This paper presents quasi-circular initial data for 729 
different configurations (mass-ratio, separation, spins) 
of binary black holes. This study covers for the first 
time a full 7-dimcnsional voluminous region of parame- 
ter space of generic non-eccentric BBH inspirals, rather 
than just lower-dimensional subspaces like aligned spin 
binaries (e.g. [13 )• The mass ratio varies between 1 and 
8 for these simulations, 620 binaries are spinning, most 
of these with generic, precessing spins of dimensionless 
magnitude up to 0.5. 

The orbital eccentricity of all 729 configurations is it- 
eratively reduced using techniques developed in earlier 
work [m . l66l . [67j . To remove or ease the computational 
burden of iterative eccentricity removal, we introduce 
fitting formulas that predict low-eccentricity orbital pa- 
rameters. Fitting the 2-dimensional non-spinning sec- 
tor spanning mass-ratio and spin requires 10 fitting pa- 
rameters each in the formulae for orbital frequency and 
radial velocity. Extending this fit to the 8-dimensional 
spinning-sector (mass-ratio, separation, two spin vectors) 
required 224 fitting parameters each. Perhaps surpris- 
ingly, despite the comparatively low spins considered here 
X ^ 0.5, higher order corrections to the spin incorporated 
via Eqs. p5gp - p5np noticeably improve the quality of 
the fit, cf. Fig. [HI We also provide a fitting formula 
for the expected number of orbits during the inspiral of 
low-eccentricity initial data, again for generic precessing 
spins and mass-ratios. 

The fitting formulas allow to achieve BH-BH config- 
urations of desired initial separation (or desired initial 



orbital frequency or desired number of orbits) with an 
eccentricity of ~ 10~^ without having to perform any it- 
erative runs. If lower eccentricity is desired, it is always 
possible to refine by iterative eccentricity removal. 

Based on these fits, we develop a technique to pre- 
dict the eccentricity of initial-data parameters without 
having to evolve the initial data at all. This technique 
is based on the deviation of the initial-data parameters 
from those for low-eccentricity (the latter determined ei- 
ther from our bank of 729 configurations, or computed 
from our low-eccentricity fits), cf. Eqs. (|^T|) - P5)) . It al- 
lows us to estimate the orbital eccentricity that would 
have resulted from post-Newtonian initial data parame- 
ters without constructing initial data for those configu- 
rations. The results are shown in Figs. [5] and [TU] demon- 
strating that our fitting formulas yield orbital eccentric- 
ity about two orders of magnitude smaller than post- 
Newtonian expansions. 

The ability to estimate eccentricity without an evo- 
lution is also useful in another scenario: To obtain an 
evolution with desired eccentricity, one can simply per- 
turb the low-eccentricity fits and use Eqs. (|4T|) - (|43| to 
determine the eccentricity of the perturbed initial-data. 
This procedure allows to construct initial data of desired 
eccentricity without having to perform any intermediate 
evolutions. This procedure is demonstrated in Sec. |Vl 

During the construction of BBH initial-data, free pa- 
rameters like initial orbital frequency and initial radial 
velocity must be chosen. Early in the evolution, the 
space-time relaxes to a quasi-stationary steady state, 
changing the black hole par amet ers slightly and emitting 
a pulse of "junk radiation" [10l[. This relaxation causes 



orbital frequency and radial velocity (and also black hole 
masses and spins, albeit to a smaller degree H^]) to de- 
viate from the corresponding parameters specified in the 
initial-data. These drifts depend on the precise type of 
initial data evolved (e.g. conformally flat conformal-thin 
sandwich data, as here; super-posed Kcrr-Schild data; 
puncture data) and may exhibit a dependence on black 
hole parameters different from the usual post-Newtonian 
terms. These drifts in parameters, induced by specific 
choices of the initial-data construction, may very well be 
the reason why our fitting formulas require terms that are 
not present in post-Newtonian expansions, cf. Eqs. (j3ip 
and ([5^ . Because these drifts — as well as the coordinate 
systems used when constructing initial data — are differ- 
ent for each class of initial-data, we expect that a similar 
fitting effort is necessary for superposed Kerr-Schild data, 
which allows access to higher black- hole spins. 

Out of the 729 configurations, we evolve 101 simula- 
tions fully through inspiral. Only the number of orbits 
of each inspiral was used here in order to write fitting 
formulas for any spinning configurations. Because of the 
wealth of information in these simulations, they will form 
the basis of a large number of future investigations into 
periastron-advance, black hole remnant properties, ana- 
lytic template modelling, and gravitational wave data- 
analysis efforts. 
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Appendix: Tables 



TABLE I: Parameters of the runs evolved in Sq. The first column is the label for each configuration (we 
refer to individual runs as "Sc21"), the following columns show the mass ratio q, the separation r/m, the initial 
orbital frequency O^, the dimensionless expansion factor d, the spin components of the first hole xa and the spin 
of the second hole xb, the eccentricity eo, the final orbital frequency flf, the number of orbits Nf between the 
initial frequency Qi and the final frequency il/. The last column denotes the number of orbits between the initial 
orbital frequency until the orbital frequency 0.05/m. 
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1.4463 


-3.60 


2.4 


12. 


19.94 


18.50 


7 


1.5 


16 




















1.4427 


-3.64 


13.0 


13. 


21.1 


19.50 


8 


1.5 


18 




















1.2199 


-2.15 


4.5 


13. 


28.9 


27.34 


9 


1.5 


14 


0.5 

















1.7367 


-5.14 


0.9 


6.3 


15.5 


14.65 


10 


1.5 


16 


0.5 


0.833 














1.4481 


-5.11 


4.3 


7.5391 


18.8 


18.00 


11 


1.5 


16 


0.5 


0.167 














1.4384 


-2.94 


0.4 


7.6682 


22.4 


21.16 


12 


1.5 


16 


0.5 


0.5 














1.4427 


-3.64 


0.4 


9.8 


20.9 


- 


13 


1.5 


16 


0.5 


1 














1.4490 


-4.76 


0.6 


7.5 


18.5 


17.74 


14 


1.5 


16 


0.5 

















1.4379 


-2.72 


1.2 


7.1 


22.5 


21.38 


15 


1.5 


17 


0.5 


1 














1.3286 


-3.82 


0.7 


7.7 


22.0 


21.27 


16 


1.5 


18 


0.5 


0.5 














1.2192 


-2.11 


3.3 


9.4 


28.9 


25.65 


17 


1.5 


19 


0.5 


1 














1.1321 


-2.26 


4.3 


7.6 


30.3 


26.08 


18 


1.5 


16 


0.5 


0.5 





0.5 


0.5 


1 


1.4424 


-3.31 


0.8 


5.6 


20.0 


19.60 


19 


1.5 


16 


0.5 


0.5 





0.5 


0.5 





1.4414 


-3.34 


4.6 


5.5 


20.1 


19.72 


20 


1.5 


16 


0.5 


0.5 





0.5 


0.5 


-0.5 


1.4534 


-3.60 


2.7 


6.8 


20.5 


19.65 


21 


1.5 


16 


0.5 


0.5 





0.5 


0.5 


0.5 


1.4533 


-3.61 


1.5 


7.8 


20.7 


19.64 


22 


1.5 


16 


0.5 


0.5 





0.5 


1 





1.4455 


-5.02 


2.8 


5.7 


18.9 


18.54 


23 


1.5 


16 


0.5 


0.5 





0.5 








1.4388 


-2.62 


2.7 


6.6 


21.6 


20.70 


24 


1.5 


16 


0.5 


1 





0.5 








1.4446 


-3.48 


1.2 


7.3 


19.7 


18.84 


25 


1.5 


16 


0.5 








0.5 


0.5 


1 


1.4370 


-2.74 


2.2 


6.8 


22.4 


21.40 


26 


1.5 


16 


0.5 








0.5 


0.5 





1.4370 


-2.74 


2.2 


6.8 


22.4 


21.40 


27 


1.5 


16 


0.5 








0.5 


0.5 


-0.5 


1.4485 


-2.57 


1.4 


5.5 


21.8 


21.40 


28 


1.5 


16 


0.5 








0.5 


0.5 


0.5 


1.4485 


-2.57 


1.4 


5.5 


21.8 


21.40 


29 


1.5 


16 


0.5 








0.5 


1 





1.4402 


-3.93 


0.2 


6.8 


21.2 


20.32 


30 


1.5 


16 


0.5 








0.5 








1.4344 


-1.56 


0.6 


7.1 


23.72 


22.38 


31 


1.5 


16 


0.320 


0.667 


0.142 


0.150 


0.732 


0.148 


1.4463 


-3.60 


2.2 


11. 


20.2 


18.80 


32 


3 


14 




















1.7427 


-4.69 


21.2 


8.3 


17.6 


16.36 


33 


3 


12 


0.5 


0.833 














2.1817 


-11.5 


4.4 


8.3 


9.3 


8.36 


34 


3 


12 


0.5 


0.167 














2.1522 


-7.72 


0.4 


9.2 


13.7 


11.71 


35 


3 


12 


0.5 


0.5 














2.1644 


-8.89 


1.6 


7.9 


11.4 


10.09 


36 


3 


12 


0.5 


1 














2.1846 


-12.3 


9.4 


10. 


9.1 


8.08 


37 


3 


12 


0.5 

















2.1508 


-7.74 


0.6 


9.3 


14.1 


11.95 


38 


3 


14 


0.5 


0.833 














1.7538 


-5.91 


0.6 


9.6 


15.3 


14.24 


39 


3 


14 


0.5 


0.167 














1.7348 


-4.24 


1.1 


8.8 


20.6 


18.69 


40 


3 


14 


0.5 


0.5 














1.7427 


-4.69 


2.1 


9.6 


18.1 


16.54 



Continued 



were performed on the GPC supercomputer at the SciNet 
HPC Consortium. SciNet is funded by: the Canada 
Foundation for Innovation under the auspices of Com- 
pute Canada; the Government of Ontario; Ontario Re- 
search Fund - Research Excellence; and the University of 
Toronto. Results obtained in this paper were produced 
using the Spectral Einstein Code (SpEC) 



TABLE I - Continued 



ScN 


q 


r/m 


XA 


Oa/tt 


4>a/t^ 


XB 








lO-'d 


10* en 




Nf 


A^O.05 


41 


3 


14 


0.5 


1 














1.7559 


-6.15 


0.3 


9.9 


14.9 


13.88 


42 


3 


14 


0.5 

















1.7338 


-3.95 


0.5 


8.7 


21.0 


18.99 


43 


3 


16 


0.5 


0.5 














1.4420 


-2.83 


3.4 


10. 


26.1 


24.63 


44 


3 


16 


0.5 


1 














1.4494 


-4.70 


29.8 


10. 


22.1 


21.13 


45 


3 


17 


0.5 


0.5 














1.3227 


-2.43 


3.5 


8.5 


30.7 


29.33 


46 


3 


17 


0.5 

















1.3175 


-2.00 


3.1 


6.1 


33.8 


32.64 


47 


3 


18 


0.5 


1 














1.2251 


-2.06 


5.4 


9.1 


31.4 


24.98 


48 


3 


14 


0.5 


0.5 





0.5 


0.5 


1 


1.7423 


-4.18 


1.7 


2.7 


11.4 


- 


49 


3 


14 


0.5 


0.5 





0.5 


0.5 





1.7409 


-4.27 


3.0 


2.7 


11.4 


- 


50 


3 


14 


0.5 


0.5 





0.5 


1 





1.7445 


-5.78 


4.2 


2.8 


11.7 


- 


51 


3 


14 


0.5 


0.5 





0.5 








1.7391 


-3.23 


4.4 


2.6 


11.1 


- 


52 


3 


14 


0.5 


1 





0.5 


1 





1.7583 


-7.93 


2.7 


5.9 


13.6 


13.23 


53 


3 


14 


0.5 








0.5 


0.5 


1 


1.7325 


-3.86 


1.7 


5.6 


19.6 


18.99 


54 


3 


14 


0.5 








0.5 


0.5 





1.7324 


-3.73 


1.8 


5.6 


19.6 


18.98 


55 


3 


14 


0.5 








0.5 


1 





1.7346 


-6.27 


6.7 


7.7 


19.9 


18.34 


56 


3 


14 


0.5 








0.5 








1.7302 


-2.86 


4.9 


5.3 


19.9 


19.56 


57 


5 


12 




















2.1686 


-5.96 


37.7 


9.1 


14.9 


13.32 


58 


5 


14 




















1.7433 


-3.50 


2.5 


9.7 


23.4 


21.69 


59 


5 


15 




















1.5809 


-2.84 


4.9 


8.3 


28.3 


26.78 


60 


5 


9.5 


0.323 


0.663 


-0.176 











3.0132 


-19.9 


48.2 


5.9 


5.6 


4.99 


61 


5 


9.5 


0.483 


0.642 


-0.212 











3.0132 


-19.9 


0.1 


8.5 


4.7 


4.73 


62 


5 


9.5 


0.5 


0.644 


-0.213 











3.0132 


-19.9 


9.4 


4.5 


4.2 


4.26 


63 


5 


12 


0.5 


0.644 


-0.213 











2.1755 


-9.02 


3.3 


2.0 


7.8 


7.75 


64 


5 


15 


0.5 


0.5 














1.5787 


-2.87 


21.6 


2.1 


13.5 




65 


5 


15 


0.5 


1 














1.5906 


-4.38 


34.2 


9.4 


22.9 


21.97 


66 


5 


15 


0.5 

















1.5695 


1.62 


42.0 


5.6 


31.7 


30.93 


67 


8 


13 




















1.9345 


-2.98 


2.9 


6.8 


24.9 


23.56 


68 


8 


13 


0.5 


1 














1.9536 


-4.21 


5.1 


7.0 


18.7 


17.99 


69 


8 


13 


0.5 

















1.9221 


-4.76 


37.8 


6.4 


31.7 


29.41 



TABLE IL Parameters of the runs evolved in S^. The first column labels the binary configuration, 
q = mA/niB denotes the mass-ratio, r/m, Qi and a denote initial separation, orbital frequency and expansion 
factor, respectively, xa and xb are the dimensionless spin components of the first and the second hole, and en 
is the orbital eccentricity. The final orbital frequency is given by and Nj denotes the number of orbits Nf 
between the initial frequency fii and the final frequency / . The last column denote the number of orbits between 
the initial orbital frequency until the orbital frequency 0.05/m. 



ScN 


1 


r/m 


XA 




(^a/t^ 


XB 


e 


b/tt 


0s /tt 


lO^mr^i 


10 




lO^en 


WmQ.j 


Nf 


A^O.05 


70 


1 


07 


15 





205 


0.965 


0.091 





385 





229 


-0.185 


1 


5821 


-5 


98 





8 


9.5 


17 


21 


15 


83 


71 


1 


08 


15 





106 


0.400 


0.590 





212 





340 


-0.191 


1 


5817 


-5 


40 





7 


13. 


17 


73 


16 


00 


72 


1 


08 


15 





154 


0.352 


-0.364 





470 





646 


-0.960 


1 


5838 


-4 


56 





2 


6.5 


16 


09 


15 


37 


73 


1 


10 


15 





494 


0.818 


-0.028 





142 





152 


-0.364 


1 


5854 


-5 


85 





5 


6.0 


15 


40 


14 


88 


74 


1 


12 


15 





269 


0.083 


0.507 





076 





474 


-0.061 


1 


5795 


-4 


35 





5 


12. 


18 


18 


16 


41 


75 


1 


12 


15 





377 


0.109 


-0.712 





458 





725 


-0.935 


1 


5809 


-3 


45 





5 


8.1 


17 


18 


15 


97 


76 


1 


12 


15 





344 


0.028 


0.492 





352 





240 


-0.931 


1 


5761 


-2 


79 





5 


11. 


19 


12 


17 


19 


77 


1 


12 


15 





243 


0.517 


-0.496 





269 





625 


-0.792 


1 


5866 


-4 


26 


2 





11. 


16 


91 


15 


43 


78 


1 


17 


15 





481 


0.958 


0.542 





370 





943 


0.220 


1 


5919 


-7 


31 


2 


3 


8.6 


14 


44 


13 


57 


79 


1 


26 


15 





267 


0.415 


-0.156 





311 





128 


-0.386 


1 


5797 


-4 


97 


1 





11. 


18 


42 


16 


66 


80 


1 


27 


15 





073 


0.786 


0.362 





463 





452 


0.106 


1 


5829 


-1 


75 





5 


7.5 


16 


88 


15 


84 


81 


1 


29 


15 





357 


0.424 


-0.907 





452 





345 


0.677 


1 


5839 


-6 


09 


2 


9 


7.8 


17 


82 


16 


59 


82 


1 


34 


15 





079 


0.568 


-0.523 





200 





853 


-0.820 


1 


5844 


-5 


45 





3 


11. 


16 


95 


15 


53 


83 


1 


35 


15 





140 


0.094 


-0.122 





413 





091 


0.430 


1 


5777 


-3 


31 





2 


10. 


18 


91 


17 


18 


84 


1 


36 


15 





089 


0.005 


-0.975 





356 





568 


-0.404 


1 


5866 


-6 


51 


2 


9 


11. 


17 


69 


16 


09 


85 


1 


38 


15 





169 


0.531 


-0.097 





455 





235 


-0.152 


1 


5799 


-5 


88 


1 





7.9 


17 


89 


16 


67 


86 


1 


55 


15 





082 


0.820 


-0.171 





205 





825 


-0.052 


1 


5847 


-5 


74 





4 


11. 


17 


20 


15 


76 


87 


1 


63 


15 





114 


0.521 


-0.665 





325 





101 


-0.587 


1 


5806 


-3 


49 


1 


2 


9.8 


18 


57 


17 


05 


88 


1 


63 


15 





295 


0.516 


0.635 





179 





060 


0.310 


1 


5832 


-5 


15 


1 


7 


12. 


18 


50 


16 


78 


89 


1 


63 


15 





171 


0.786 


-0.032 





155 





544 


-0.983 


1 


5843 


-5 


15 





2 


11. 


17 


46 


15 


99 


90 


1 


64 


15 





115 


0.435 


0.808 





343 





730 


-0.361 


1 


5861 


-6 


82 


1 


4 


11. 


17 


65 


16 


17 


91 


1 


68 


15 





365 


0.148 


-0.283 





217 





269 


-0.911 


1 


5772 


-3 


53 





8 


12. 


20 


04 


18 


03 
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TABLE II - Continued 



ScN 


q 


r/m 








XB 


93 /-K 


4>b/tt 




lO^'d 


lO-^en 


lO^mfi/ 


Nj 


iVo.05 


92 


1.70 


15 


0.429 


0.470 


0.487 


0.112 


0.958 


-0.874 


1.5880 


-5.20 


4.9 


12. 


18.22 


16.60 


93 


1.76 


15 


0.406 


0.530 


-0.810 


0.252 


0.764 


0.476 


1.5866 


-2.88 


1.5 


11. 


17.79 


16.33 


94 


1.82 


15 


0.403 


0.326 


-0.511 


0.146 


0.590 


-0.953 


1.5828 


-4.13 


2.4 


11. 


19.44 


17.67 


95 


1.82 


15 


0.438 


0.658 


-0.446 


0.124 


0.253 


0.690 


1.5891 


-5.74 


2.9 


9.5 


17.56 


16.29 


96 


1.84 


15 


0.208 


0.364 


0.875 


0.160 


0.420 


0.640 


1.5815 


-4.99 


0.9 


11. 


19.08 


17.40 


97 


1.91 


15 


0.375 


0.545 


-0.854 


0.110 


0.012 


-0.107 


1.5826 


-2.71 


0.6 


11. 


18.76 


17.17 


98 


1.92 


15 


0.451 


0.082 


-0.887 


0.443 


0.563 


-0.652 


1.5834 


-H3.02 


4.6 


8.6 


20.34 


18.73 


99 


1.94 


15 


0.312 


0.873 


0.292 


0.222 


0.080 


-0.934 


1.5853 


-4.26 


0.8 


11. 


17.88 


16.46 


100 


1.96 


15 


0.057 


0.810 


0.274 


0.090 


0.408 


-0.113 


1.5827 


-4.68 


0.4 


12. 


18.76 


17.14 


101 


1.97 


15 


0.190 


0.564 


-0.917 


0.372 


0.911 


0.727 


1.5851 


-5.78 


0.5 


9.5 


17.82 


16.56 



TABLE III: Parameters ai^j of pNS- TABLE IV: Parameters of kns- 



i 


j = 


j = l 


i = 2 


J =3 




i 


i = o 


J = l 


j = 2 


j = 3 





-2.23773 


89.0053 


-907.744 


3194.70 







0.999768 


-2.63274 


4.08714 


-2.67330 


1 


16.3786 


-321.623 


1589.74 







1 


-0.0126939 


0.343292 


1.75169 





2 


-25.8886 


285.810 










2 


0.0445636 


-0.975111 








3 


7.94034 













3 


0.0157333 












TABLE V: Parameters 



of Shis. 



J = 1 


i = 2 


J = 3 


j = 4 


5316.71 


-66338.2 


275375. 


-380348. 


32067.7 


-376132. 


1.48805x10" 


-1.93024x10" 


14274.6 


-177525. 


736095. 


-1.01733x10" 


-1843.34 


9864.44 


1329.56 


-4958.02 


109.978 


517.973 


-7645.17 


15290.1 


23169.4 


-275526. 


1.08011x10" 


-1.42663x10" 


7819.45 


-101310. 


434622. 


-617796. 


-42885.3 


532586. 


-2.20338x10" 


3.03600x10" 


1027.50 


-27767.7 


103643. 


-218964. 


-109193. 


1.29526x10'^ 


-5.05535x10" 


6.62736x10" 


-55234.7 


704761. 


-2.98370x10" 


4.19269x10" 


-92749.8 


1.02247x10'^ 


-4.02065x10" 


5.29669x10" 


-7921.47 


97674.7 


-405463. 


565779. 


-103301. 


1.22607x10'* 


-4.82157x10" 


6.34451x10" 


-82960.2 


1.03434x10** 


-4.29063x10" 


5.92217x10" 


-496416. 


5.85007x10'^ 


-2.32014x10'^ 


3.02622x10' 


-217655. 


2.71089x10" 


-1.12563x10'' 


1.55773x10' 


25018.3 


-393524. 


2.05655x10" 


-4.12184x10" 


28030.5 


-368855. 


1.57451x10" 


-2.19102x10" 


-412124. 


4.87348x10" 


-1.90352x10'' 


2.49538x10' 


-128584. 


1.65610x10" 


-7.06852x10" 


1.00034x10' 


661241. 


-8.21574x10" 


3.40040x10'' 


-4.68715x10' 


-24548.5 


478963. 


-1.79194x10" 


3.37544x10" 


1.69467x10'' 


-2.01047 xlO'^ 


7.86851x10'' 


-1.03283x10* 


917140. 


-1. 16445 xlO'^ 


4.90819x10' 


-6.86985x10' 


-289050. 


3.65366x10" 


-1.08819x10' 


8.95028x10" 


46482.3 


-583937. 


2.51586x10" 


-3.69108x10" 


1.74966x10'^ 


-2.07098x10'^ 


8.12790x10' 


-1.06573x10* 


424444. 


-5.28782x10" 


2.19181x10' 


-3.02304x10' 


2.54951 xlO'^ 


-3.01226 xlO'^ 


1.19551x10* 


-1.56397x10* 


1.08832x10'^ 


-1.35750 xlO'^ 


5.64446x10' 


-7.82131x10' 


538865. 


-4.65423x10" 


1.08665x10' 


-650921. 


-134766. 


1.75054x10" 


-7.36970x10" 


1.00998x10^ 


2.38473x10" 


-2.81363 xlO'^ 


1.09819x10* 


-1.43508x10* 


690728. 


-8.84984x10" 


3.76020x10' 


-5.30045x10' 


-3.34460 xlO'^ 


4.15737x10^ 


-1.72135x10* 


2.37354x10* 


149314. 


-2.61396x10" 


1.01824x10' 


-1.80892x10' 


-8.70071 xlO'^ 


1.03298x10* 


-4.05413x10* 


5.32989x10* 



0,1 
(0,2) 
0,3 

0,4; 

(0,5) 
0,6 

0,7; 

(0,8) 
(0,9) 
(0,10) 
(0,11) 
0,12 
0,13 
(0,14) 

(1,1) 
(1,2) 
(1,3) 
(1,4) 
(1,5) 
(1,6) 
(1,7) 
(1,8) 
(1,9) 
(1,10) 

(1,11) 
(1,12) 
(1,13) 
(1,14) 
(2,1) 
(2,2) 
(2,3) 
(2,4) 
(2,5) 
(2,6) 
(2,7) 
(2,8) 
(2,9) 
12M 
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.ed 


{i, k) 


i = 1 


i = 2 


J = 3 


j = 4 


(2,11) 


-4.98684x10" 


6.30344x10'' 


-2.64631x10** 


3.69064x10* 


(2,12) 


362328. 


2.89435x10" 


-6.26526x10'' 


1.54943x10* 


(2,13) 


-296859. 


3.76776x10" 


-1.62996x10'^ 


2.38893x10' 


(2,14) 


-9.73227x10'' 


1.15044x10** 


-4.51208x10* 


5.90608x10* 


(3,1) 


-709184. 


8.82898x10" 


-3.65705x10' 


5.04053x10' 


(3,2) 


-4.35305x10'' 


5.15001x10^ 


-2.04343x10* 


2.67729x10* 


(3,3) 


-1.78215x10'' 


2.22596x10^ 


-9.26725x10' 


1.28564x10* 


(3,4) 


908964. 


-1.60025x10^ 


8.51391x10'' 


-1.47147x10* 


(3,5) 


62120.2 


-864159. 


3.64740x10" 


-4.74819x10" 


(3,6) 


-4.48307x10" 


5.28627x10^ 


-2.06464x10* 


2.69548x10* 


(3,7) 


-1.20683x10" 


1.53949x10^ 


-6.51614x10' 


9.15436x10' 


(3,8) 


5.53612x10" 


-6.88414 xlO'^ 


2.85135x10* 


-3.93290x10* 


(3,9) 


-265327. 


4.52185x10" 


-1.85463x10' 


3.22321x10' 


(3,10) 


1.48479x10^ 


-1.76472x10* 


6.94463x10* 


-9.14535x10* 


(3,11) 


8.83120x10" 


-1.11213x10* 


4.65325x10* 


-6.46984x10* 


(3,12) 


-6.70535x10" 


5.32904x10^ 


-6.84880x10' 


-1.07762x10* 


(3,13) 


1.00833x10" 


-1.27036 xlO'^ 


5.39248x10' 


-7.69297x10' 


(3,14) 


1.77463x10^ 


-2.09671x10* 


8.22363x10* 


-1.07571x10^ 



TABLE VI: Parameters di^j^u of 5ps. 



(i, k) 


j = 1 


j = 2 


J = 3 


i = 4 


(0,1) 


-45382.3 


302204. 


-153169. 


-1.30349x10" 


(0,2) 


-253943. 


1.90739x10" 


632048. 


-7.50294x10" 


(0,3) 


-1.31299x10" 


1.54603x10' 


-6.03927x10' 


7.82728x10' 


(0,4) 


-1.20248x10" 


8.98755x10" 


-2.09367x10' 


2.51332x10' 


(0,5) 


550699. 


-6.20436x10" 


2.42440x10' 


-3.28133x10' 


(0,6) 


587207. 


-7.79345x10" 


3.28891x10' 


-4.78587x10' 


(0,7) 


-234908. 


2.08401x10" 


-5.31451x10" 


2.96473x10" 


(0,8) 


529948. 


-4.20533x10" 


7.57153x10" 


2.97058x10" 


(0,9) 


677458. 


-1.23094x10' 


5.30880x10' 


-9.37766x10' 


(0,10) 


4.56142x10" 


-5.15215x10' 


2.05792x10* 


-2.57620x10* 


(0,11) 


-5.37955x10" 


8.86716x10' 


-4.55622x10* 


7.47340x10* 


(0,12) 


-1.23435x10' 


1.01134x10* 


-3.15628x10* 


3.51274x10* 


(0,13) 


-1.45409x10" 


1.69206x10' 


-6.86296x10' 


9.64437x10' 


(0,14) 


5.15330x10" 


-5.79718x10' 


2.18460x10* 


-2.69339x10* 


(1,1) 


-45382.3 


302204. 


-153169. 


-1.30349x10" 


(1,2) 


-253943. 


1.90739x10" 


632048. 


-7.50294x10" 


(1,3) 


-1.31299x10" 


1.54603x10' 


-6.03927x10' 


7.82728x10' 


(1,4) 


-1.20248x10" 


8.98755x10" 


-2.09367x10' 


2.51332x10' 


(1,5) 


550699. 


-6.20436x10" 


2.42440x10' 


-3.28133x10' 


(1,6) 


587207. 


-7.79345x10" 


3.28891x10' 


-4.78587x10' 


(1,7) 


-234908. 


2.08401x10" 


-5.31451x10" 


2.96473x10" 


(1,8) 


529948. 


-4.20533x10" 


7.57153x10" 


2.97058x10" 


(1,9) 


677458. 


-1.23094x10' 


5.30880x10' 


-9.37766x10' 


(1,10) 


4.56142x10" 


-5.15215x10' 


2.05792x10* 


-2.57620x10* 


(1,11) 


-5.37955x10" 


8.86716x10' 


-4.55622x10* 


7.47340x10* 


(1,12) 


-1.23435x10' 


1.01134x10* 


-3.15628x10* 


3.51274x10* 


(1,13) 


-1.45409x10" 


1.69206x10' 


-6.86296x10' 


9.64437x10' 


(1,14) 


5.15330x10" 


-5.79718x10' 


2.18460x10* 


-2.69339x10* 


(2,1) 


-45382.3 


302204. 


-153169. 


-1.30349x10" 


(2,2) 


-253943. 


1.90739x10" 


632048. 


-7.50294x10" 


(2,3) 


-1.31299x10" 


1.54603x10' 


-6.03927x10' 


7.82728x10' 


(2,4) 


-1.20248x10" 


8.98755x10" 


-2.09367x10' 


2.51332x10' 


(2,5) 


550699. 


-6.20436x10" 


2.42440x10' 


-3.28133x10' 


(2,6) 


587207. 


-7.79345x10" 


3.28891x10' 


-4.78587x10' 


(2,7) 


-234908. 


2.08401x10" 


-5.31451x10" 


2.96473x10" 


(2,8) 


529948. 


-4.20533x10" 


7.57153x10" 


2.97058x10" 


(2,9) 


677458. 


-1.23094x10' 


5.30880x10' 


-9.37766x10' 


(2,10) 


4.56142x10" 


-5.15215x10' 


2.05792x10* 


-2.57620x10* 


(2,11) 


-5.37955x10" 


8.86716x10' 


-4.55622x10* 


7.47340x10* 


(2,12) 


-1.23435x10' 


1.01134x10* 


-3.15628x10* 


3.51274x10* 


(2,13) 


-1.45409x10" 


1.69206x10' 


-6.86296x10' 


9.64437x10' 


(2,14) 


5.15330x10" 


-5.79718x10' 


2.18460x10* 


-2.69339x10* 
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TART 



K VT - CnnUn 



imL. 



.7 = 3 



(3,1) 
(3,2) 
(3,3) 
(3,4) 
(3,5) 
(3,6) 
(3,7) 
(3,8) 
(3,9) 
(3,10) 
(3,11) 
(3,12) 
(3,13) 

(3M 



-45382.3 

-253943. 
-1.31299x10® 
-1.20248x10® 

550699. 

587207. 

-234908. 

529948. 

677458. 
4.56142x10® 
-5.37955x10® 
-1.23435x10^ 
-1.45409x10® 
5.15330x10® 



302204. 
1.90739x10® 
1.54603x10^ 
8.98755x10® 
-6.20436x10® 
-7.79345x10® 
2.08401x10® 
-4.20533x10® 
-1.23094x10^ 
-5.15215x10^^ 
8.86716x10^ 
1.01134x10* 
1.69206x10^ 
-5.79718x10'' 



-153169. 

632048. 
-6.03927x10'' 
-2.09367x10' 
2.42440x10'' 
3.28891x10' 
-5.31451x10® 
7.57153x10® 
5.30880x10' 
2.05792x10* 
-4.55622x10* 
-3.15628x10* 
-6.86296x10' 
2.18460x10* 



7 



-1.30349x10" 
-7.50294x10® 
7.82728x10' 
2.51332x10' 
-3.28133x10' 
-4.78587x10 
2.96473x10® 
2.97058x10® 
-9.37766x10' 
-2.57620x10* 
7.47340x10* 
3.51274x10* 
9.64437x10' 
-2.69339x10* 



Advanced ligo, http://lhocds.ligo- 

wa.caltech.edu:8000/advligo/GWINC. 

L. scientific collaboration, The science of Isc research, [19] 
URL http : //www. ligo . org/ science/ overview. php 
J. Abadie et al. (The LIGO Scientific Collaboration [20 
and the Virgo Collaboration, the Virgo Collaboration) [21 
(2012), 1201.5999. 

F. Acernese, M. Alshourbagy, P. Amico, F. Antonucci, [22 
S. Aoudia, et al.. Class. Quant. Grav. 25, 184001 (2008). 
K. Somiya and the KAGRA Collaboration, Class. 
Quantum Grav. 29, 124007 (2012). [23 
F. Pretorius, Class. Quantum Grav. 22, 425 (2005). 
J. Centrella, J. G. Baker, B. J. Kelly, and J. R. van 
Meter, Rev. Mod. Phys. 82, 3069 (2010). [24 
H. P. Pfeiffer, Numerical simulations of compact object 
binaries (2012), 1203.5166. 

A. Buonanno, Y. Pan, J. G. Baker, J. Centrella, B. J. [25 

Kelly, et al., Phys. Rev. D 76, 104049 (2007), 0706.3732. 

P. Ajith, S. Babak, Y. Chen, M. Hewitson, B. Krishnan, [26' 

et al., Phys. Rev. D 77, 104017 (2008), 0710.2335. 

T. Damour and A. Nagar, Phys. Rev. D 79, 081503 [27 

(2009), 0902.0136. 

A. Buonanno, Y. Pan, H. P. Pfeiffer, M. A. Scheel, [28 
L. T. Buchman, et al., Phys. Rev. D 79, 124028 (2009), 
0902.0790. [29 
P. Ajith, M. Hannam, S. Husa, Y. Chen, B. Bruegmann, 

N. Dorband, D. Mueller, F. Ohme, D. PoUney, C. Reis- [30 
swig, et al. (2009), 0909.2867. 

Y. Pan, A. Buonanno, L. T. Buchman, T. Chu, [31 
L. E. Kidder, et al., Phys.Rev. D81, 084041 (2010), 
0912.3466. 

B. Aylott et al.. Class. Quantum Grav. 26, 165008 [32 
(2009). 

L. Santamaria, F. Ohme, P. Ajith, B. Briigmann, 
N. Dorband, M. Hannam, S. Husa, P. Mosta, D. PoUney, 

C. Reisswig, et al., Phys. Rev. D 82, 064016 (2010). [33 
P. Ajith, M. Boyle, D. A. Brown, B. Brug- 
mann, L. T. Buchman, et al.. Class. [34 
Quantum Grav. 29, 124001 (2012), URL 
[http : //stacks . iop . org/0264-9381/29/i=12/a=124001 . 
B. Aylott, J. G. Baker, W. D. Boggs, M. Boyle, P. R. [35 



Brady, et al., Class.Quant.Grav. 26, 165008 (2009), 
0901.4399. 

P. C. Peters and J. Mathews, Phys. Rev. 131, 435 
(1963). 

P. C. Peters, Phys. Rev. 136, B1224 (1964). 

F. Ohme, Class.Quant.Grav. 29, 124002 (2012), 
1111.3737. 

J. G. Baker, S. T. McWilliams, J. R. van Meter, J. Cen- 
trella, D.-I. Choi, B. J. Kelly, and M. Koppitz, Phys. 
Rev. D 75, 124024 (2007). 

B. Briigmann, J. A. Gonzalez, M. Hannam, S. Husa, 
U. Sperhake, and W. Tichy, Phys. Rev. D 77, 024027 

(2008) . 

M. Campanelli, C. Lousto, P. Marronetti, and 
Y. Zlochower, Phys.Rev. Lett. 96, 111101 (2006), gr- 
qc/0511048. 

M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 
Rev. D 74, 041501 (R) (2006), gr-qc/0604012. 
S. Husa, J. A. Gonzalez, M. Hannam, B. Briigmann, and 
U. Sperhake, Class. Quantum Grav. 25, 105006 (2008). 
L. Boyle, M. Kesden, and S. Nissanke, Phys. Rev. Lett. 
100, 151101 (2008), arXiv:0709.0299 [gr-qc]. 
J. G. Baker, M. Campanelli, F. Pretorius, and Y. Zlo- 
chower, Class. Quantum Grav. 24, S25 (2007). 
T. Bode, D. Shoemaker, F. Herrmann, and I. Hinder, 
Phys. Rev. D 77, 44027 (2008). 

T. Chu, H. P. Pfeiffer, and M. A. Scheel, Phys. Rev. D 
80, 124051 (2009), 0909.1313. 

M. A. Scheel, M. Boyle, T. Chu, L. E. Kidder, K. D. 
Matthews and H. P. Pfeiffer, Phys. Rev. D 79, 024003 

(2009) , arXiv:gr-qc/0810.1767. 

M. Hannam, S. Husa, J. G. Baker, M. Boyle, B. Brueg- 
mann, T. Chu, N. Dorband, F. Herrmann, I. Hinder, 
B. J. Kelly, et al., Phys. Rev. D 79, 084025 (2009), 
arXiv:0901.2437. 

G. Lovelace, M. Scheel, and B. Szilagyi, Phys. Rev. D 
83, 024010 (2011), 1010.2777. 

M. Hannam, S. Husa, F. Ohme, D. Miiller, and 
B. Briigmann, Phys. Rev. D 82, 124008 (2010), 
arXiv: 1007.4789. 

G. Lovelace, M. Boyle, M. A. Scheel, and 



18 



B. Szilagyi, Class. Quant. Grav. 29, 045003 (2012), [61 
arXiv:1110.2229. 

L. T. Buchman, H. P. Pfeiffer, M. A. Scheel, and B. Szi- [62; 

lagyi. Simulations of unequal mass binary black holes 

with spectral methods (2012), 1206.3015. [63; 

M. Hannam, S. Husa, F. Ohme, and P. Ajith, Phys. 

Rev. D 82, 124052 (2010). 

P. Ajith, S. Babak, Y. Chen, M. Hewitson, B. Krishnan, [64 

A. M. Sintes, J. T. Whelan, B. Briigmann, P. Diener, 

N. Dorband, et aL, Phys. Rev. D 79, 129901 (2009). 

P. Ajith, M. Hannam, S. Husa, Y. Chen, B. Bruegmann, [65 

N. Dorband, D. Mueller, F. Ohme, D. PoUney, C. Reis- 

swig, et al., ArXiv e-prints (2009), 0909.2867. 

A. L. Tiec, A. H. Mroue, L. Barack, A. Buonanno, [66 

H. P. Pfeiffer, et al., Phys.Rev.Lett. 107, 141101 (2011), 

1106.3278. 

A. H. Mroue, L. E. Kidder, and S. A. Teukolsky, Phys. [67 
Rev. D 78, 044004 (2008). 

A. Taracchini, A. Buonanno, E. Barausse, M. Boyle, 
T. Chu, G. Lovelace, H. P. Pfeiffer, and M. A. Scheel [ 
(2012), 1202.0790. [69^ 
M. Boyle, A. Buonanno, L. E. Kidder, A. H. Mroue, 
Y. Pan, et al., Phys. Rev. D 78, 104020 (2008), [70 
0804.4184. 

Y. Pan, A. Buonanno, M. Boyle, L. T. Buchman, L. E. 
Kidder, H. P. Pfeiffer, and M. A. Scheel, Phys. Rev. D [71 
84, 124052 (2011), 1106.1021. 

M. Campanelli, C. O. Lousto, Y. Zlochower, and [72 
D. Merritt, Astrophys. J. Lett. 659, L5 (2007). [73^ 
M. Campanelli, C. O. Lousto, Y. Zlochower, B. Krish- 
nan, and D. Merritt, Phys. Rev. D 75, 064030 (2007), [74 
gr-qc/0612076. [75^ 
P. Schmidt, M. Hannam, S. Husa, and P. Ajith, Phys. 
Rev. D 84, 024046 (2011), arxiv: 1012.2879. [76; 
OShaughnessy, R. and Vaishnav, B. and Healy, J. and 
Meeks, Z. and Shoemaker, D., Efficient asymptotic [77 
frame selection for binary black hole spacetimes using 
asymptotic radiation (2011), arXiv:1109.5224. [78 
M. Boyle, R. Owen, and H. P. Pfeiffer, Phys. Rev. D [79^ 
84, 124011 (2011), arXiv:1110.2965. 

R. Sturani, S. Fischetti, L. Cadonati, G. Guidi, J. Healy, [80 

et al. (2010), 1012.5172. [81 

R. Sturani, S. Fischetti, L. Cadonati, G. Guidi, J. Healy, [82 

et al., J.Phys.Conf.Ser. 243, 012007 (2010), 1005.0551. 

R. O'Shaughnessy, J. Healy, L. London, Z. Meeks, [83 

and D. Shoemaker, Phys.Rev. D85, 084003 (2012), [84 

1201.2113. 

A. Buonanno, G. B. Cook, and F. Pretorius, Phys. Rev. [85 
D 75, 124018 (2007), gr-qc/0610122. [; 
H. P. Pfeiffer, D. A. Brown, L. E. Kidder, L. Lindblom, 
G. Lovelace, and M. A. Scheel, Class. Quantum Grav. [87 
24, S59 (2007). 

J. G. Baker, J. R. van Meter, S. T. McWilliams, J. Cen- [ 
trella, and B. J. Kelly, Phys. Rev. Lett. 99, 181101 
(2007). [89 

C. Lincoln and C. Will, Phys. Rev. D 42, 1123 (1990). 
T. Damour and G. Schafer, Nuovo Cimento Soc. Ital. 

Fis. 101 B, 127 (1988). [90; 

T. Damour, A. Gopakumar, and B. R. Iyer, Phys. Rev. 

D 70, 064028 (2004). [91 

C. Konigsdorffer and A. Gopakumar, Phys. Rev. D 73, 

124012 (2006), gr-qc/0603056. [92; 

R.M. Memmesheimer, A. Gopakumar and G. Schafer, 

Phys. Rev. D 70, 104011 (2004). [93; 



E. Berti, S. Iyer, and C. M. Will, Phys. Rev. D 74, 
061503 (2006), gr-qc/0607047. 

T. Mora and C. Will, Phys. Rev. D 66, 101501(R) 

(2002) . 

S. Husa, M. Hannam, J. A. Gonzalez, U. Sperhake, 
and B. Briigmann, Phys. Rev. D 77, 044037 (2008), 
0706.0904. 

M. Campanelli, C. O. Lousto, H. Nakano, and Y. Zlo- 
chower, Phys. Rev. D 79, 84010 (2009), arXivigr- 
qc/0808.0713. 

A. H. Mroue, H. P. Pfeiffer, L. E. Kidder, and 
S. A. Teukolsky, Phys. Rev. D 82, 124016 (2010), 
arXiv: 1004.4697 [gr-qc]. 

A. Buonanno, L. E. Kidder, A. H. Mroue, H. P. Pfeif- 
fer, and A. Taracchini, Phys.Rev. D83, 104034 (2011), 
1012.1549. 

M. Boyle, D. A. Brown, L. E. Kidder, A. H. Mroue, H. P. 

Pfeiffer, M. A. Scheel, G. B. Cook, and S. A. Teukolsky, 

Phys. Rev. D 76, 124038 (2007). 

W. Tichy and P. Marronetti (2010), 1010.2936. 

M. Purrer, S. Husa, and M. Hannam, Phys.Rev. D85, 

124051 (2012), 1203.4258. 

L. T. Buchman, H. P. Pfeiffer, M. A. Scheel, and B. Szi- 
lagyi, Simulations of non- equal mass black hole binaries, 
in preparation. 

D. Brown and P. Zimmerman., Phys. Rev. D 81, 024007 

(2010), arXiv:gr-qc/0909.0066. 

J. W. York, Phys. Rev. Lett. 82, 1350 (1999). 

H. P. Pfeiffer and J. W. York, Phys. Rev. D 67, 044022 

(2003) . 

G. B. Cook, Phys. Rev. D 65, 084003 (2002). 

G. B. Cook and H. P. Pfeiffer, Phys. Rev. D 70, 104016 

(2004) . 

M. Caudill, G. B. Cook, J. D. Grigsby, and H. P. Pfeif- 
fer, Phys. Rev. D 74, 064011 (2006). 

H. P. Pfeiffer, L. E. Kidder, M. A. Scheel, and S. A. 
Teukolsky, Comput. Phys. Commun. 152, 253 (2003). 
http : //www. black-holes . org/SpEC . html 

L. Lindblom, M. A. Scheel, L. E. Kidder, R. Owen, and 

0. Rinne, Class. Quantum Grav. 23, S447 (2006). 

H. Friedrich, Commun. Math. Phys. 100, 525 (1985). 

D. Garfinkle, Phys. Rev. D 65, 044029 (2002). 

C. Gundlach, J. M. Martin-Garcia, G. Calabrese, and 

1. Hinder, Class. Quantum Grav. 22, 3767 (2005). 
O. Rinne, Class. Quantum Grav. 23, 6275 (2006). 

O. Rinne, L. Lindblom, and M. A. Scheel, Class. Quan- 
tum Grav. 24, 4053 (2007). 

J. M. Stewart, Class. Quantum Grav. 15, 2865 (1998). 
H. Friedrich and G. Nagy, Commun. Math. Phys. 201, 
619 (1999). 

J. M. Bardeen and L. T. Buchman, Phys. Rev. D 65, 
064037 (2002). 

B. Szilagyi, B. Schmidt, and J. Winicour, Phys. Rev. D 
65, 064015 (2002). 

G. Calabrese, J. PuUin, O. Reula, O. Sarbach, and 
M. Tiglio, Commun. Math. Phys. 240, 377 (2003), gr- 
qc/0209017. 

B. Szilagyi and J. Winicour, Phys. Rev. D 68, 
041501(R) (2003). 

L. E. Kidder, L. Lindblom, M. A. Scheel, L. T. Buch- 
man, and H. P. Pfeiffer, Phys. Rev. D 71, 064020 (2005). 
L. T. Buchman and O. C. A. Sarbach, Class. Quantum 
Grav. 23, 6709 (2006). 

L. T. Buchman and O. C. A. Sarbach, Class. Quantum 



19 



Grav. 24, S307 (2007). 
[94] D. Gottlieb and J. S. Hesthaven, J. Comput. Appl. 

Math. 128, 83 (2001), ISSN 0377-0427. 
[95] J. S. Hesthaven, Appl. Num. Math. 33, 23 (2000). 
[96] S. Ossokme, L. E. Kidder, and H. P. Pfeiffer, 

Precession-tracking coordinates for binary black hole 

simulations, in preparation. 
[97] B. Szilagyi, L. Lindblom, and M. A. Scheel, Phys. Rev. 



D 80, 124010 (2009), 0909.3557. 
[98] L. E. Kidder, Phys. Rev. D 52, 821 (1995). 
[99] B. Garcia, G. Lovelace, L. E. Kidder, M. Boyle, S. A. 
Teukolsky, et al. (2012), 1206.2943. 
[100] G. Lovelace, R. Owen, H. P. Pfeiffer, and T. Chu, Phys. 

Rev. D 78, 084017 (2008). 
[101] G. Lovelace, Class. Quantum Grav. 26, 114002 (2009). 



